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A two-dimensional analysis of the potential in charge coupled devices 
is presented. It is assumed that there are no mobile minority carriers, that 
the plate separation is zero, and that the plate voltage does not vary with 
time. The depletion layer approximation is used to linearize the equations, 
which are then solved exactly with the use of Fourier series. Both surface 
and buried channel devices are analyzed. These solutions can typically be 
evaluated on a computer in less than a tenth of the time it takes to obtain 
a solution by the method of finite, differences. The solutions obtained here 
provide an important tool for the designer of charge coupled devices. In 
addition to describing the method of obtaining the solutions, we evaluate 
them to show the effects of a number of different design parameters, and 
compare the cost of these solutions with the cost of obtaining finite difference 
solutions. 

I. INTRODUCTION AND SUMMARY 

The recent invention 1 - and development of charge coupled devices 
(CCD's) has led to renewed interest in the mathematical analysis of 
MIS-type structures. Ideally, one would like to solve the nonlinear 
equations describing the three-dimensional motion of charge as a func- 
tion of the time-varying plate voltages. So far no one has succeeded in 
doing this for even the simplest geometries. For the most part, one- 
dimensional static models have been solved which yield only qualita- 
tive information about the behavior of such devices. A much more 
sophisticated, one-dimensional, time-varying model of a CCD has been 
developed and analyzed by Schryer and Strain. 3 

A static, two-dimensional model of a CCD has also been studied by 
Amelio 4 using finite difference techniques. He calculated the potential 
distribution in a two-dimensional model in the absence of mobile 
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charge and with given static plate potentials. The results of this calcu- 
lation are of great interest. The use of finite difference techniques in 
these cases has drawbacks, however. In even the relatively simple 
geometries considered so far, it is expensive to obtain reasonably ac- 
curate solutions for the potentials, and for more complicated devices, 
it soon becomes prohibitively expensive. Furthermore, as we shall 
show, even for simple geometries it is difficult to obtain accurate ex- 
pressions for the fields from the finite difference solutions for the 
potentials. 

In this paper, we show that when the plates on a CCD are close 
enough together so that they can be assumed to be abutted, and when 
the depletion layer approximation can be used, 8 the resulting linearized 
model can be solved analytically. These solutions can then be evalu- 
ated cheaply and quickly on a computer. This analysis will be valid for 
both surface and buried channel CCD's with an arbitrary number of 
plates. In a separate paper, we will show that these solutions can be 
used to obtain solutions for the potential in a CCD when there are 
gaps between the plates. 6 

In Section II we write down the equations describing the model and 
put them into appropriate dimensionless form. We then introduce the 
depletion layer approximation which linearizes the equations and dis- 
cuss conditions under which this approximation is valid. 

This paper has two main purposes: to show the behavior of the po- 
tentials and fields in a CCD and to demonstrate techniques by which 
these potentials and fields can be calculated cheaply and accurately. 
In Section III we present a discussion of how the solutions depend on 
the various parameters defining the devices. 

In Section IV we derive the solution of the linearized potential equa- 
tions. The reader interested only in the physical design of CCD's can 
skip the rest of the paper. 

In Section V we discuss in some detail the solution by finite differ- 
ence methods of the exact, nonlinear equations describing a surface 
CCD. Our purpose in doing this is twofold. We wish to show the diffi- 
culties involved in obtaining an accurate solution cheaply, especially 
if an accurate knowledge of the fields is required. Secondly, we want 
accurate solutions of the exact problem to compare with the analytic 
solutions of the linearized problem. 

Finally, in Section VI we compare in detail some solutions of the 
exact problem obtained by finite differences with the corresponding 
analytic solutions of the linearized equations. It is shown that in many 
cases of interest the solutions of the linearized problem provide excel- 
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lent approximations to the true potential and much more accurate ap- 
proximations to the gradient of the potential than can be obtained from 
the finite difference solutions. Furthermore, the solutions of the linear- 
ized equations are at least an order of magnitude cheaper to obtain 
than are the finite difference solutions for any reasonable accuracy. 

II. DERIVATION OF THE EQUATIONS 

A surface CCD 1 consists of a layer of silicon covered with a thin 
insulating layer of silicon dioxide, and on top of the oxide layer, a 
sequence of closely spaced electrodes. Such a device is shown schemati- 
cally in Fig. 1 with some typical dimensions indicated. Mobile charge 
trapped at the oxide-semiconductor interface is transferred from plate 
to plate by appropriately changing the potential of the plates. We con- 
sider the case where the substrate is n-type silicon and the mobile 
charges are injected holes. In this case, the plate potentials must be 
negative. Our analysis can be modified in an obvious way to describe 
the case where the substrate is p-t3 r pe silicon and the mobile charges 
are electrons. 

Some losses are introduced by the trapping of the mobile charges by 
surface states at the oxide-semiconductor interface. Smith and Boyle 7,8 
have proposed a solution to this problem by inserting between the oxide 
and the substrate an additional layer of p-type silicon, thus forming 
a buried channel CCD. The p-layer is kept completely ionized, which 




0.1/im 



Fig. 1 — A schematic diagram of a surface CCD. 
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Fig. 2 — A schematic diagram illustrating plate overlap structure. 

causes the potential minimum to occur near the center of the p-layer. 
Thus the mobile positive charge stays at the potential minimum, safely 
away from the surface traps. 

It is desirable to have the plates as close together as possible in order 
to minimize the transit time of the mobile charge between plates. The 
minimum plate separation presently obtainable from photolithog- 
raphy is ~3-5 Mm, but a plate overlap structure, 89 as shown in Fig. 
2, or undercut isolation 10 allows for plate separation of ~0.1 /mi. 

We propose to study the static potential in either a surface or buried 
channel CCD with plate overlap structure, in the absence of mobile 
charge. We begin by noting that since the length in the ^-direction of 
each plate is much greater than its width in the ^-direction, near the 
center of the plates (z = 0) the field is essentially two-dimensional. 
Hence we will treat the problem as two-dimensional. We assume that 
the plates are zero distance apart. Since in the overlap structure there 
should be little flux leakage between the plates, we feel this is a 
reasonable approximation. We make the additional assumption that 
the bottom substrate is infinitely thick. The field can penetrate into 
the substrate little beyond a depletion depth, and since for typical 
voltages the depletion depth ranges from 7 /jm to 20 j*m, and the thick- 
ness of the substrate in a typical device is 100 tan, this is a very rea- 
sonable approximation. Finally, we assume the structure is periodic in 
the x-direction, which in the usual mode of operation is an excellent 

approximation. 

We begin by denning the boundary value problem describing a 

buried channel device. In all that follows, starred quantities have 

rationalized MKS dimensions; unstarred quantities, except for a few 

obvious physical parameters, are dimensionless. In the strip 

^ x* ^ L*, let <p\ denote the potential in the oxide layer, ^ y* 

^ hi; <p' 2 the potential in the p-layer, h\ ^ y* ^ hi; and <p' 3 the po- 
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Fig. 3 — A schematic diagram of one cell of a three-phase, buried channel CCD. 



tential in the substrate, h* 2 ^ y* < oo (see Fig. 3). Then 

V 2 ,<fil = 0, (1) 

V;^ = j 2 N- A (y*) (2) 

v>5 - - ^ (1 - exp M/kT)). (3) 

In eqs. (1) through (3), V*. is the two-dimensional Laplacian; -e is the 
charge of an electron; N A (y*) is the acceptor number density in the 
p-region; N' D is the donor number density in the n-region; ei and e 2 
are the permittivity of the oxide and silicon, respectively; k is Boltz- 
man's constant; and T is the absolute temperature. The conditions 
under which eq. (3) can be expected to be valid are discussed in Ref. 
11. 

In the usual method of fabricating a CCD, the substrate donor 
number density, N' D , is a constant, independent of position. However, 
in a buried channel CCD, the p-layer is formed by diffusing the ac- 
ceptor ions into the n-type Si, and so N* A (y*) is typically a function of 
y* — h\, the distance from the oxide surface. In many cases we have 
the representation 12 

N' D 



N' A {y*) = Caexp 



-(fcHS)'* 



-N' D 



(4) 
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where C& is the number density of acceptor ions at the upper surface of 
the Si. The average value, N* A) of N' A (y) is easily shown to be 



"2 — "I J hi 



Vtt Cs 



erfWCW) -N' D (5) 



(7) 



2 Vfri(c.;/AQ 

where erf (x) is the error function. I3 

Before writing down the boundary conditions, we introduce dimen- 
sionless variables. Define the Debye length 

X D = (e 2 kT/e 2 N' n )K (6) 

Then normalize all lengths with respect to \p, 

x = x*/\ D , y = y*/\ D , L = L*/X D , 

h k = h' k /X D , (k = 1, 2), 
and define 

<Pa{x, y) = e&x*, y*)/kT, (a = 1, 2, 3), (8) 

*{y) = N A (y*)/N' D , v = 6i/«. (9) 

Equations (1) through (3) become 

VVi = 0, (10) 

vvi = <y), ( n ) 

Wa = e« - 1. (12) 

The boundary conditions can be written now as follows : For ^ x ^ L, 

n (s,0) = F(.r), (13) 

*(*, fa) = i^s, *i)i * ^ («. *0 = ^ (*i *0 + Q(«), (14) 

M*, h>) = *> s (*, U, ^ fe W = ^ («, W, (is) 

«(*,•) -0, ( 16 ) 

and 

<?(0, ?/) = ¥>(£, 2/) 

^ (0, y) = d -£ (L, y), 0g ? ,<*. (17) 

In (13), 7(x) is a given, periodic function, assuming on each electrode 
the constant voltage of the electrode; and in (14), Q(x) is a known, 



CCD POTENTIAL 



675 



periodic surface charge density, which may include deliberately im- 
planted charges. 14 

For future use, we record the expression for a(y) when Na(v*) is 
given by (4). If 

C s = C;/N' D , (18) 

then 

a(y) = Cs exp j - ( jjLZ^i p»C, s l - 1, (19) 

Vir C s erf (^l(nCs) 



and 



jr-nCt 



- 1. 



(20) 



The equations for the surface CCD are essentially the same, except 
that the p-layer is eliminated. In what follows, we will only give the 
analysis for the buried channel CCD. The results for the surface CCD 
can be obtained from those for the buried channel CCD by setting 
it = 0, Ai = ho, and <^ 3 = *pz. 

We now introduce the important depletion layer approximation. 5 
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Fig. 4 — A schematic diagram of the depletion layer approximation for one cell 
of a three-phase, buried channel CCD. 
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In most cases of interest, p 3 (x, h 2 ) « — 1 for ^ x ^ L. For example, 
in a typical buried channel CCD, <p 3 (x, h 2 ) ~ - 1000 [and in a typical 

surface CCD, <p 2 (x, hi) > 150]. Thus, for y — h 2 small and positive, 

e *3 ^ o. However, for y » h 2 , \ cp 3 \ « 1, and e* 3 = 1 + <p3. There is 
thus some curve y = R(x) such that for h 2 ^ y ^ R(x),e™ — 1 » — 1. 
The region h 2 ^ y ^ 7?(.r) is the depletion region. For R(x) < y, we 
have e v " 3 — 1 ~ tp*. If R{x) varies but little about its average value, 
R, these remarks suggest that we replace eqs. (10) through (12) by the 
system of linear equations 



VVi (x, y) = 0, 
V 2 fa (x, y) = <r(y), 
V^ 3 (x,y) = -1, 
VVi (x, V) = 1A4 (a:, y), 



< y < k h (21) 

hr<y < h 2 , (22) 

/* 2 < y < h 3 = h 2 + i?, (23) 

/i 3 = ht + i? < y < oo (24) 



where ^ 3 is the potential in h 2 ^ y ^ h 3 and ^ 4 is the potential in 
hs ^ y < °° . (See Fig. 4.) In addition to ^ h \J/ 2 , and \J/ 3 satisfying 
boundary conditions (13) through (15), we have the boundary condi- 
tions, for ^ x ^ L, 



fa (x, h 3 ) = fa (x, h 3 ), -r- 3 (x, h 3 ) = -^~ (x, h 3 ), 



d^i 



Bfa 



By 



By 



(25) 
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Fig. 5 — The potential <p*(x*, y*) plotted as a function of x* in a surface CCD for 
y* = 0.2 (im, 2.275 pm, and 4.35 pm. The 45-/im plates are alternately at -5 V and 
— 10 V, and the oxide thickness is hf = 0.2 ^m. 
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Fig. 6— The field - (d<p*/dx*) (x* y*) plotted as a function of x* in a surface CCD 
for y* = 0.2 M m. The 45-^m plates are alternately at —5 V and - 10 V, and the oxide 
thickness ia h* = 0.2 M m. The dashed curve is a plot of —(d<p*/dx*) for the same 
device obtained from a finite difference calculation. 



iM*, °°) = 0, 



(20) 



and the 0„ (a = 1,2, 3, 4) all satisfy (17). 

It has been shown that for a one-dimensional version of this problem, 
the choice 5 

R = - (l + h 2 - ft, + h -l ) + [(/, 2 - fc, + ^ - 1 - 2V 

) d i ^ hl + ^)^ &di ] i (2?) 






yields a solution which approximates the solution of the nonlinear 
problem very accurately in the region h x g y <J /; 2 . Furthermore, the 
solution in this region is quite insensitive to the choice of R. Since an 
accurate knowledge of the potential is only necessary in the p-layer for 
the buried channel CCD and near the oxide-semiconductor interface 
for the surface CCD, we feel this approximation is well justified. In 
this two-dimensional problem, we determine R from (27) by letting V 
and Qss be the averages of V(x) and Q(x): 



Vo = jj'' V(x)dx, Qss = IJ^ Q{x)dx. 



(28) 
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Fig. 7 — The potential <p*(x*, y*) plotted as a function of x* in a surface CCD for 
y* = 0.2 fim, 2.275 /xm, 4.35 ftm, and 6.425 /*m. The 5-/zm plates are alternately at 
-5 V and —10 V, and the oxide thickness is h? = 0.2 ^m. 

In Section VI we will present a comparison of the solution of the 
linearized equations with the solution of the nonlinear equations ob- 
tained by finite difference methods. This will confirm for this example 
that the approximate solutions are accurate as claimed. 

III. GENERAL BEHAVIOR OF THE POTENTIALS AND FIELDS 

In this section we present graphical representations of the potentials 
and fields for both surface and buried channel CCD's for a number of 
design parameters. The graphs were obtained by evaluating the ana- 
lytic expressions for the solutions, derived in Section IV, of eqs. (13) 
through (15), (17), and (21) through (26). 

In all cases, we assume that the doping in the n-type substrate is 
N D - 10 14 /cm 3 , that e 2 /e„ = 12, where e is the permittivity of free 
space, that ei/« g = i and that Q(x) = 0, i.e., there is no trapped or 
implanted charge at the semiconductor-oxide interface. Then at 
T = 300°K, the value of \ D defined in (6) is X D = 4.15 X lO" 5 cm. 
Also, in all the examples presented here, we have used the factor (kT/e) 
= 0.025 V to convert dimensionless potentials to volts, and the factor 
(kT/e\ D ) = 600 V/cm to convert dimensionless fields to V/cm. 

We consider first the effect of plate width in surface devices. The first 
pair of graphs illustrate a surface CCD with 45-mui plates, the second 
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pair a surface CCD with 5-/xm plates, and the third pair a surface CCD 
with 1.5-fxm plates. The oxide layer in each of these CCD's is 0.2 ^ni 
thick, and the voltages on the plates are alternately —5 V and — 10 V. 
These examples show the storage mode, and as a result there is no 
asymmetry to introduce a perferred direction of flow for the holes. 
However, they do exhibit the effects of plate width, and are easy to 
compare with finite difference calculations. In Figs. 5, 7, and 9, <p* is 
plotted along the oxide-semiconductor interface (y* = h\) and along 
the lines y* = h\ + 5\ D , h\ + \0\ D , and h\ + lb\ D inside the sub- 
strate. In Figs. (5, 8, and 10, -(d<p*/dx*) = E' x is plotted along the 
oxide-semiconductor interface. The dashed curve is the field calculated 
by finite difference methods. The discrepancy between the two curves 
will be discussed in Section V. 

In all three cases, the peak field available for moving positive charge 
from the left-hand plate to the center plate is about 4.8 X 10 4 V/cm. 
However, in the 45-jum plate device, the field penetrates only about 
7 /jiii under the plate from the edge, leaving most of the region under 
the plate field free. This would clearly be a very poor CCD. On the 
other hand, in the 1.5-fim device, there are substantial fields under the 
whole plate. These graphs show that if field penetration under the 
plates were the sole criterion, the narrower the plates the better. How- 
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Fig. 8— The field — (d^*/d.r*)(.r*, //*) plotted as a function of x* in a surface CCD 
for //* = 0.2 fim. The 5-/im plates are alternately at —5 V and — 10 V, and the oxide 
thickness is /if = 0.2 /um. The dashed curve is a plot of —(dip*/dx*) for the same 
device obtained from a finite difference calculation. 
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Fig. 9— The potential v*(x*, y*) plotted as a function of x* in a surface CCD for 
y* = 0.2 fim, 2.275 pm. 4.35 ^m, and 6.425 »m. The 1.5-/xm plates are alternately at 
—5 V and —10 V, and the oxide thickness is h{ = 0.2 /im. 

ever, recent work by Tompsett 16 has shown that in surface CCD's 
the difference in losses of ones and fat zeros due to surface states be- 
comes greater as the plate width decreases. His work shows that this 
puts a lower bound on plate widths in the neighborhood of 5 jum. How- 
ever, Fig. 8 shows that for 5-jum plates there is still considerable field 
penetration under the plates. 

Our calculations show that increasing (decreasing) the thickness of 
the oxide layer decreases (increases) the peak values of the fields, but 
does not materially affect the penetration of the fields under the plates. 

We next consider buried channel CCD's. As for surface devices in 
general, the narrower the plates the better as far as field penetration is 
concerned. However, the plate width is apt to be determined by cur- 
rent photolithography tolerances, so this is a parameter not easily 
varied. In addition, if the plates are too narrow, the charge-carrying 
capacity of the CCD becomes very small. 

Instead of considering the effects of plate width, we examine what 
happens for a given plate width if the thickness of the p-type layer is 
varied. We consider first a three-phase, buried channel CCD with 
5-vm plates. The plates are at -5 V, - 10 V, and - 15 V, so charge is 
to be moved from under the - 10 V plate to under the - 15- V plate. 
The thickness of the oxide layer is h\ = 0.1 nm. The doping profile 
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Fig. 10— The field - (d<p*/dx*) (»*, ;/*) plotted as a function of x* in a surface CCD 
for y* = 0.2 /*m. The 1.5-jum plates are alternately at —5 V and — 10 V, and the oxide 
thickness is /t,* = 0.2 pm. The dashed curve is a plot of — {d<p*/dx*) for the same 
device obtained from a finite difference calculation. 



in the p-type layer is assumed to be given by eq. (4) with 
Cs = 4.6 X 10 15 /cm 3 (Cs = 46). This corresponds to an average 
number density of acceptor atoms of 2 X 10 15 /cm 3 . The remaining 
physical parameters are as described at the beginning of the section. 
In Fig. 11 we plot the channel field, E* x = —(dtp*/dx*), (that is, the 
field at the potential minimum in the p-region) as a function of x* for 
three different p-type layer thicknesses: h' 2 — h\ = 0.1 jum, 2 jim, and 
4 urn. In Fig. 12 we plot the corresponding channel potentials, <p*, (that 
is, the value of the potential at the potential minimum in the p-layer). 
The CCD with h' 2 — h\ = 0.1 /xm is essentially a surface device. As the 
thickness of the p-layer is increased, the minimum value of the field 
under the center of the — 10-V plate increases at first, while the peak 
value of the field decreases. Eventually, as the thickness of the p-layer 
is increased, the channel will be so far below the plates that the channel 
fields will start decreasing to zero. Thus, in terms of field penetration, 
there appears to be an optimal p-layer thickness. From Fig. 12, it is 
clear that as the p-layer gets thicker, the channel potential curve 
flattens out, and so the charge-carrying capacity of the CCD decreases. 
We have also studied the effects of varying the doping in the p-layer 
(i.e., varying Cs). The behavior of the fields is relatively insensitive to 
changes in C s . 
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Fig 1 1 — The channel field — (d <p*/di*) plotted as a function of .r* along the chan- 
nel for three buried channel CCD's. The 5-miii plates are at -5 V, - 10 V, and - 15 V, 
and the oxide thickness is /i,* = 0.1 urn. C% = 4.6 X 10 15 /cm 3 , and the p-layers are 
0.1 /xm, 2 ixm, and 4 nm thick, respectively. 

An operating buried channel CCD has been reported 16 in which the 
gaps between the plates have been filled with a resistive material so 
that the potential drop between the plates is essentially linear. The 
butted plate model can be easily adapted to describe this. In Figs. 13 
and 14 we show the channel fields and potentials in such a CCD. This 
is a three-phase CCD with 10-/zm plates and 5-Mm gaps. The voltages 
on the plates are -5 V, - 10 V, and - 15 V and in the gaps the voltage 
varies linearly from one plate to the next. The thickness of the oxide 
layer is h[ = 0.1 /xm. The doping profile in the p-type layer is given by 
(4) with C's = 4.6 X 10 15 /cm 3 (C s = 46). The remaining physical 
parameters are as described at the beginning of this section. In Fig. 
13 we plot the channel field, El = -(dip*/dx*), as a function of x* 
for two different p-layer thicknesses, h' 2 - h\ = 3 pm and 7 urn. In 
Fig. 14 we plot the corresponding channel potentials, <p*. Again, it is 
seen that the p-layer thickness is a sensitive parameter in terms of field 
penetration and charge-carrying capacity, and there is undoubtedly 
an optimal thickness. In Figs. 15 and 16 we plot the same quantities 
for a three-phase CCD which is identical to the one of Figs. 13 and 14, 
except that the gap spacing is zero. The devices seem to have essentially 
the same fields and charge-carrying capacities. 
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IV. ANALYTIC SOLUTION OF THE LINEARIZED EQUATIONS 

In this section we briefly derive the solution of eqs. (21) through (24) 
subject to the boundary conditions (13) through (15), (17), (25), and 
(26). We assume as given the Fourier series expansions of V(x) and 
Q(x): 



V(x) = 2«o + Z (fl« cos X„ x + b„ sin \„ x), 

n-l 



(29) 

(30) 
(31) 



Q(x) = U o + Z (f » eos X„ x + £„ sin X„ x) 

where 

\„ - (2nw)/L. 

Since in most cases of interest V(x) and Q(x) are either piece wise con- 
stant or linear functions, it is trivial to obtain the coefficients of these 
series. 

Since ^(.r, ij) must be periodic in x with period L, we can expand the 
solution in each of the four regions in a series of the form 

\P(x, y) = \A {y)+ f (A „ (y) cos A„ x + B„ (y) sin X„ x). (32) 
»— i 

On substituting expressions of the form (32) for \p into (21) through (24) 




Fig. 12 — The channel potential <p* plotted as a function of x* along the channel 
for three buried channel CCU's. The 5-/im plates are at —5 V, —10 V, and — 15 V, 
and the oxide thickness is /i, = 0.1 nm. C% = 4.6 X 10 15 /cm 3 , and the p-layers are 
are 0.1 pm, 2 pm, and 4 pm thick, respectively. 
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Fig. 13 — The channel field - (d<p*/dx*) plotted as a function of x* along the chan- 
nel for two buried channel CCD's. The 10-^m plates are at —5 V, - 10 V, and - 15 V 
and are separated by 5-Mm gaps in which the potential varies linearly between plates. 
kf = 0.1 Mm, C% = 4.6 X KF/cm', and the p-layers are 3 Mm and 7 pm thick, 
respectively. 
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Fig. 14 — The channel potential <p* plotted as a function of x* along the channel for 
two buried channel CCD's. The 10-mhi plates are at -5 V, - 10 V, and -15 V and 
are separated by 5-^m gaps in which the potential varies linearly between plates. 
hf = 0.1 Mm, C% = 4.6 X 10 16 /cm 3 , and the p-layers are 3 Mm and 7 Mm thick, 
respectively. 
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Fig. 15 — The channel field — {d^/dx*) plotted as a function of x* along the chan- 
nel for two buried channel CCD's. The 10-/*m plates are at —5 V, — 10 V, and — 15 V, 
h* = 0.1 ftra, C's = 4.6 X 10'Vcm 3 , and the p-layers are 3 /an and 7 /an thick, 
respectively. 

and equating to zero the coefficients of cos X„ x and sin \ n x,n = 0, 1, 2, 
• • • , we obtain an uncoupled system of second-order, constant- 
coefficient, ordinary differential equations from which the A„ (y) and 
B„ (y) can be determined simply. Each A n (y) and B n (y) is the sum of 
two linearly independent solutions and thus each involves two con- 
stants of integration which must be determined by making use of the 
boundary conditions (13), (14), (15), (17), (25), and (26). Since the 
Fourier series representing the solutions must be equal term by term 
at the boundaries, this yields a simple set of linear algebraic equations 
for the unknown constants of integration. These equations can be 
solved explicitly, yielding the integration constants as linear functions 
of the coefficients a„ and b„, and £ „ and £„, of the Fourier series for V(x) 
and Q{x) given in (29) and (30). The algebra involved is elementary 
but involved, and we only record the final answer here. 



Let 



F» (x) = a n cos X„ x + b„ sin \ n x, 



*n (z) = f n cos A,, x + £„ sin X„ x, 
so that we can write 



(33) 
(34) 



V(x) = \a + Z Fn Or), 

n— I 



(35) 
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Fig. 16 — The channel potential <p* plotted as a function of x* along the channel for 
two buried channel CCD's. The 10-^m plates are at -5 V, -10 V. and -15 V, 
/if =0.1 nm, C% = 4.6 X lO'Vcm 3 , and the p-layers are 3 nm and 7 pm thick, 
respectively. 



Q(X) =hto+ Z *»(*)■ 
n-1 



(36) 



Furthermore, let 

tf ± = 1 d= I/17, (37) 

At = 1 ± (1 + K 2 )*, (38) 

M n {y) = \E + A+ + E_A-e- n »V' 3 - hi) \e- x '"> 

+ {E-A+e- 2X "' 11 + E + A-e- n »'' 3 \e*"", (39) 

L„ (?/) = 2\A+e- Xn,J + A-e-^' 2/ '3-^}- (40) 

Then we can write 
ft (x, y) = (Ad + B) + (Cfi„ + D)2/ 

+ f If (x) m " iy) + * ( c ) L " (/ * l)sinhx " y l (4i) 

+ ky J ' ( j M„ (0) + " W M„ (0) ^ I ' < ' 
*2 (x, y) = \ \a„ (1 + h 3 - h 2 ) - (h* - h,) 2 

- (a + 2h 2 - 2h)(y - h 2 ) + 2 J" (y - £V(£)^1 
sinh A„ hi 1 L„ (y) 



+ Z \Fn(x) + *„(*) 



nX, 



M„ (0) ' 



(42) 



CCI) POTENTIAL 687 

* (•>•>.</) - K&o (l + An - y) ~(y- h a )*l 

+ f (,.„ + ., „-*£*)£$, (43) 

+ 41 {*.(,) + #.<,)!«**•*! 1 

ex p e- VTTxfa - *i) - X- fti] ,^x 
if. (0) l j 

where 

4 =1(1 + /,,-/^ + ^), (45) 

B = J* U ~ h + j ) *(£)# - g [(Ai - /^)(//3 + h - 2fti) 



Z) = 

and 



f e + 2(fc 8 - //•,) -2[ l '<r(&dl 

J h 



+ ^ (f . + 2/ i3 - 

>7 


- 2/i 2 ) 


, (46) 


C = -1/(2,), 




(47) 


r^cf)^ 1/(211), 




(48) 



a = (fa, - B)/ii. (49) 



Equations for drfr/dx and d\f//dy can be obtained by differentiating 
eqs. (41) through (44) term by term. To obtain the equations for a sur- 
face CCD, drop equation (42) and set a = 0, h 2 = hi, and relabel ha, 
\f/ 3 , and ^4 as /?2, ^2, and ^3 in the remaining equations. 

V. FINITE DIFFERENCE SOLUTION OF THE EXACT EQUATIONS 

In this section we describe the finite difference (FD) solution 17 of 
a surface CCD described by equations (10) through (17). We will 
assume that Q(x) = 0, a(y) = 0, and hi = h 2 = h. 

The infinite region ^ x ^ L, h ^ y < 00 is replaced by the finite 
rectangle ^ x ^ L, h ^ y ^ H, with H ~2> h, and the boundary con- 
dition (16) is replaced for ^ x ^ L by 

Vi (.r, H) = 0. (50) 

This may be done, because the solution tends to zero rapidly as y ~ ► «> . 
In fact, a one-dimensional analysis 3 shows that <^2 tends to zero ex- 
ponentially in y. 

A uniform FD net is now placed over region 1 (the oxide layer) and 
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region 2 (the silicon layer) with N points in the ^-direction and M i and 
M2 points in the ^-direction in regions 1 and 2, respectively. Let 

h z = L/(N - 1), h lv = h/(M t - 1), 
ho y = (H - h)/(M 2 - 1), 

and then, for (1 £ i ^ -V), (1 ^ J ^ M«), and (a = 1,2), define 

¥>«.,-./= *.((»'- 1)A„ 0'- DW (52) 

where <?„ foy) is the exact solution of (10) through (17), and define 
^a, i,7 as the FD solution which approximates <p „,.,>. 

The FD equations are obtained as follows. The boundary condition 
(50) is replaced by the N equations 

02.i.*, = O, (lS*£JV)i (53) 

while the boundary condition (13) is replaced by the N equations 

Mil = V«i - l)h x ) = V i} (l£ig N). (54) 

If x = iio — l)h x is the edge of a plate, the jump discontinuity in 
V(x) there is handled by setting 

V io = HV io+ i + V ie „i). (55) 

At each interior mesh point V 2 v„,,,y is approximated by the standard 
five-point difference formula 17 

Vft<p a .i,j = (<Pa,i+l,j + <Pa.i-l,j ~ %0a,i,j)/hl 

+ ($a,i,3+\ + 0a,i,j-l — 2<P a ,i,i)/hly. (56) 

Equation (56) can also be used to evaluate Vl<p a ,\,j and Vl0 a ,N.j, 
(2 ^ j ^ M a — 1) by making use of the periodicity relations 

0a,-l,j = 0a,N-l.j, 0a,N+\,j = <P*,\,j- (57) 

Thus, eqs. (10) and (12) are replaced by (Mi + M 2 - 4)iV equations 

V^i.v.y = 0, (Ut^ N), (2 ^ j ^ M, - 1), (58) 

V^,,,-., = exp (0 2 ,ij) - 1, (lli^ AT), 

(2 ^ J ^ Mi - 1). (59) 

There remain the interface conditions (14). The first of these is re- 
placed by the N equations 

0i. i ,M l = 0*.i.i, (l£i£N). (60) 

To obtain an equivalent set of equations for the second condition, we 
could replace the derivatives by the first differences from each side. 
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However, it is well known that this approximation is not very good. 
This is easily seen from the prototype equation V 2 <^ = 0, where, for 
simplicity, we take h z = h„ — h. Then as is well known, 17 

<p i+lJ + ^_ lti + ?,,,•+, + <^,, y _ 1 - 4<f tJ = A*V"*« + 0(/i 4 ), (61) 
while 

*.*., - p.-.y = h ^ (a-,, y y ) + 0(A»). (62) 

Thus, the errors in the FD equations, after scaling the left-hand sides 
to have coefficients of order 1 , are out of balance. The interface condi- 
tion is l//i 2 less accurately modeled than the differential equation. This 
leads to the following scheme which balances the errors equally. We 
want d<pi/dy using only values of tp h similarly for dif-i/dy, and the 
approximation must be good to 0(h 3 ). This may be done using the 
values of ipx.ij, (Mi — 3 ^ ./' ^ Mi). Simply use the derivative of the 
cubic interpolation polynomial through these values. It is easily seen 
that 

^(Xi,h) = [llipi.i.Mt ~ 18 



<fl,i,Ml-l 

+ 9^,,,. Ul _o - 2*i, I -.j, 1 _ 1 }/(dfti,) + 0(A?»). (63) 



Similarly, 

^ (x it h) = \2vt,t, A - 9*«.<., + 18^.,-. 2 

-ll*» Mi i)/(<tt*) + 0(Al). (64) 

Then the second boundary condition (14) is replaced by the .V equa- 
tions 

1l{W01,i,Mi ~ I8<£l. ,..!/,-! + 9< 1 3l,i,.U 1 -2 — 2#l,i, Mi-*\ /h\ y 
= {202.1.4 — 9^2,i. :i + 1802, i.2 — H02,i.l\/hi v , 

(1 £ t £ N). (65) 

Equations (53), (54), (58), (59), (60), and (65) comprise {Mi + M 2 )N 
equations in the (Mi + M 2 )N unknowns #«,*,/, (a = 1, 2), (1 £ { ^ N), 
and (1 ^ j ^ M a ). From the standard FD theory, the solution of this 
set of transcendental equations differs from the solution of the true 
boundary value problem by a factor of order 0(/j 2 ). We used a nonlinear 
overtaxation scheme developed in Refs. 18 through 20 to solve the 
FD equations for 0% and standard overrelaxation methods 21 to solve 
the 0i FD equations. 

An initial estimate of the solution, 0®l Jt was obtained by computing 
the one-dimensional matching solutions as functions of y along the 
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lines x = .r„ (1 ^ i ^ N), using the methods of Ref. 5. This provides 
a solution which is a fairly good estimate under the middle of any plate 
and a very bad one near the edge of any plate. These one-dimensional 
solutions also give a good estimate of the greatest depth of the transi- 
tion region, call it y mttX . Since <ps -> exponentially for y > y max , we 
chose H = 7/, uax + 20. 

This estimated solution &% is now refined iteratively by the method 
of successive overrelaxation (SOR). The (n + l)st iterate is obtained 
from the nth as follows. For all n = 0, 1, 2, • • • , set 

0[% = V it 0i tt U = 0, (l£i£ N). (66) 

For (2 ^ i ^ N - 1) and (2 ^ j ^ M, - 1), let pfcj" be denned by 

(2K 2 + mnti£* = (*©HJ + *fc±fc)k" 9 

+ Mh+i + dSPWAS 1 . (67) 

Then set 

0$}" = CO! pftj" + (1 - udiQj ( 68 ) 

where wj is an overrelaxation parameter satisfying 1 ^ cm < 2. For 
(2 ^ « ^ JV — 1) and (2 ^ j ^ M 2 — 1), let £ 2 .,-.,• be the solution of 

(2hr 2 + 2V)0£tf" + exp (d& n ) = 1 + (#©hj + rtStVj)*^ 

Equation (69) has the form 

Aw + a* = S (70) 

where A and £ are known and A > 0. Given any approximate solution 
w° of (70), Newton's method 22 yields the sequence 

w *+i = [^ + C '"*]->[# + («,* - !)«"*! (A: = 0, 1, 2, • • •), (71) 

which converges to the solution of (69). The convergence of this scheme 
is global and quadratic because the function Aw + e"\ for A > 0, is 
a monotone increasing, convex function of iv. After solving for p|*f "» 
we set 

$&' = " 2 0ft)* + (1 - «.)*& (72) 

where 1 ^ to* < 2. 

The interface values $i, i ,M 1 and 2 ,.\i, (2 ^ « i JV — 1) are re- 
laxed by combining (60) and (65), defining 

tttP, = dSt" = W + *i l )-»Ci*5 l (180iSffi-i - WWU 

+ 2fftffi-,) + ^'(18^ - 90$, + 2^ 4 )]/(ll), (73) 
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and setting 

0W> = 0&' = <*0t&\ + (1 - »)*Rm (74) 

where 

« = |(wj + wa). (75) 

Equations (OG) through (75) describe the manner in which the in- 
terior nodes of the FD mesh are relaxed. The nodes at i = 1 and N 
involve periodicity and require more detailed study. By using the 
periodicity relations (57), equations (66) through (75) can be extended 
to the nodes at i = 1, IV. We could next do one of two things. First, we 
could treat a ,ij &nd a ,N,i as separate quantities and relax each of 
them separately using the periodicity relations (57). Then, each time 
a , ij was relaxed, this new value could be substituted for a ,n,j to 
preserve periodicity, and vice versa. This was tried but gave very poor 
convergence rates. The problem is that in treating a ,\,i an ^ 0a, nj 
as separate quantities, the same quantity (that is a .u and a ,N.i) 
gets relaxed twice rather than once in each SOR sweep. This can be 
avoided by letting a ,yj = a ,\.ii (« = 1, 2), (2 ^ j S M* — 1) and 
then relaxing only the quantities a ,ij. This produces quite acceptable 
convergence rates. 

The overrelaxation parameters a>i and o> 2 were set equal to the opti- 
mum values of these parameters for the Laplace equation on regions 
one and two respectively. These values for «i and w 2 were estimated as 
follows. Let <p a n) denote the vector of values of the nth iterate of the 
solution to V\0 a = in region a, (a = 1, 2). Define the nth residual 
vector as R a n) = $ a " + ]) — 0™. Then, starting with any initial guess 
a O) 5* 0, standard theory shows 23 that lim ||Ri , " H> ll/l|Ri ,,) ll = v« exists 
and 

w a .o P , = 2/jl + \1 - Va ] (76) 

where ||Ri n, iJ denotes the norm of the vector R^"' and is called the resid- 
ual. In practice, we calculated the residuals for n large enough so that 
?7a was obtained to the desired accuracy using the L-> norm. 

A further important point is solving the transcendental equation 
(69), which must be done at each SOR step. In the (n + l)th SOR 
sweep, the initial estimate for 02jj u in the Newton iteration (71) was 
02.u- There is no reason to compute the ^o" j"" very accurately when 
0$$j is far away from its final value. Conversely, the quadratic con- 
vergence of Newton's method means that when the error in 0j$j is 
small, one Newton iteration will produce a very good approximation 
for <?2"J 1) - For this reason, only one Newton iteration was used in solv- 
ing (69) during each SOR sweep. 
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Some theory has been developed to show that the SOR scheme we 
have outlined above converges to the true solution in some mildly non- 
linear Dirichlet problems. 1819 However, to the best of our knowledge, 
no theoretical analysis exists of the boundary value problem of this 
paper. Nevertheless, as we will demonstrate by numerical examples in 
the next section, the scheme works in practice. 

We conclude this section with a few remarks on the estimated ac- 
curacy of the FD solution. Since little is known of the general theory 
of a complicated nonlinear boundary value problem such as we are 
considering, we argue by analogy with the Dirichlet problem for 
Laplace's equation on a square. Let <p(z, y) be the true solution of 
V 2 <p = in < x, y < L, with <p specified on the boundary, let 
h = L/(N - 1), <p {lj = <p((i - l)h, (./ - l)h), and let U be the 
solution of the corresponding FD equations. Then it is well known 24 
that under reasonable conditions on the boundary values, 

\\<P - 0\U = SU P l*« "~ 00"l = 0(^ 2 )- ( 77 ) 

This relationship assumes that we know the FD solution exactly. How- 
ever, we don't know exactly, all we know are the various iterates <p {n) 
which have been calculated and the residuals ||<£<'' +1) - 0< B) |L. Now 
it is known, though not as widely as it should be, when calculating 
by the method of SOR that 25 

||0 - (B+1) IU = C(o>)||0<» +1) - M \U, (78) 

and if cot is the optimal choice of w, 

CM = 0(N). (79) 

It can also be shown, if the optimal value of u is underestimated by ten 
or fifteen percent, that C(«) = 0(iV 2 ). This means that to obtain an 
approximate solution accurate to Q(h 2 ) = 0(N~ 2 ) by the method of 
SOR, we must iterate at least until the residuals are 0(JV- 3 ), and since 
in the nonlinear problem we can only crudely estimate the optimal w, 
we should really iterate until the residuals are 0(N~*). 

As we will show by example in the next section, it is necessary to 
calculate the potentials with great accuracy if one wishes to obtain the 
fields from them with any accuracy at all by differencing them. From 
the previous paragraph, however, we have shown that this is expensive 
in even a moderately complicated problem, since then the residuals 
must be made so small. To estimate the cost of increasing the size of the 
problem or decreasing the mesh size (both equivalent to increasing N), 
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we note that typically in SOR, 26 

||0<»+I) _ 0(„)|| m « g - n |0(,V-l)| ( 80 ) 

for optimal w, while for nonoptimal w, 

||0(n+l) _ 0(»>|| M « e -»IO(JV-I)| ( gl ) 

If wc wish to specify that ||$ — <£ ( " +1) ||« ^ e, then, for optimal u>, it 
is easy to show from (78) through (80) that the number of iterations 
must be 

n = 0(N tn{N /<■)), (82) 

while for nonoptimal w, it follows from (78), (79), and (81) that 

n = 0(NHn(N/e)). (83) 

Thus, if we decrease the mesh size by one half, then from (82) or (83), 
we must double or quadruple the number of iterations to obtain the 
same accuracy. Since there are now 4 times as many mesh points, the 
time required to obtain a solution goes up by a factor of 8 or 16, de- 
pending on the knowledge of o>. 

VI. COMMENTS ON THE ACCURACY OF THE SOLUTIONS 

In Section IV we derived the Fourier series solution (FSS) of the 
linearized problem, and in Section V we outlined the finite difference 
solution (FDS) of the nonlinear problem. In this section we compare 
several of these solutions with regard to accuracy and cost. All calcu- 
lations discussed were performed on a Honeywell 6070 computer, and 
all programs were written in Fortran IV. 

We solved the nonlinear equations (10) through (17) by the method 
of finite differences for the two-phase surface CCD's discussed in Sec- 
tion III, some of whose properties are presented graphically in Figs. 
5 through 10. The three CCD's have plate widths of 45 /im, 5 Mm, 
and 1.5 pm, respectively. The plate voltages are —5 V and —10 V, 
Q(x) = 0, N D = 10 14 /cm 3 , ci/eo = 4, e 2 /e = 12, and \ D = 0.415 /im. 
In all cases, a FD net was chosen with N = 25, il/i =» 25, and Alz = 41. 

In the case of the 45-/zm-plate device, this corresponds in the dimen- 
sionless units to h x = 9, hi u = 0.02, h 2y = 1. After 273 SOR iterations, 
the residual was ~4 X 10~ 7 , the running time was 195 seconds, and 
35 K of memory was used. This should ensure that the difference be- 
tween the true FDS and the iterated solution will never exceed 
~25 X 65 X 4 X 10~ 7 = 6.5 X lO" 4 . We have calculated the FSS at 
the same mesh points, this took 33 seconds to run, and 33 K of memory 
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was used. Let tp FD (x, y) denote the FD solution of the nonlinear prob- 
lem and \f/(x, y) denote the FS solution of the linearized problem, and 

e(x, y) = | (^fd (x ,y) - $(x, y))/<p F D(x, y) \ . (84) 

Then we found that along the oxide-semiconductor interface, e(x, hi) 
< 1.14 X 10 -2 , ^ x ^ L, and five Debye lengths below this inter- 
face, e(x, ^ + 5) < 2.89 X 10~ 2 , ^ x ^ L 

We have for the 5-Mm-plate device h x = 1, h\ y = 0.02, h 2u = 1. 
After 288 SOR iterations, the residual was ~4 X 10" 7 , the running 
time was 190 seconds, and 35 K of memory was used. This should again 
ensure that the difference between the true FDS and the iterated 
solution does not exceed 6.5 X 10 -4 . For this case, we found that 
e(x, h^ < 2.8 X 10- 3 , e(x, fn + 5) < 1.4 X 10~ 3 , e(x, hi + 10) < 3.14 
X 10- 3 , e(x, In + 15) < 9.81 + 10~ 3 for ^ x ^ L. The running 
time to evaluate the FSS was 10 seconds and 33 K of memory was used. 

Finally, for the 1.5-Mm-plate device, we have h x = 0.3, h iu = 0.02, 
and h 2u = 1. After 300 SOR iterations the residual was ~6 X 10~ 4 , 
the running time was 196 seconds, and 35 K of memory was used. Note 
that in this case the residual is three orders of magnitude greater than 
in the other two cases. We found that e(x, hi) < 1.8 X 10~ 2 , the run- 
ning time to calculate the FSS was 4 seconds, and 33 K of memory was 
used. 

In Ref. 5 it was noted that as long as |^ 4 (h 3 ) + 1 1 < 10, one could 
expect the solution of the linearized problem to be a good approxima- 
tion to the solution of the nonlinear problem, at least in the p-region 
for buried channel devices or near the oxide-semiconductor interface 
for surface devices. In the examples considered here, for the 45-^m- 
plate case, -5.15 < ^ 4 (x, h 3 ) < 3.00, for the 5-Mm-plate device, 
-1.20 < ti{x, h 3 ) < -0.96, and for the 1. 5-Mm-plate case, ^ 4 (x, h 3 ) 
= —1.077. This again suggests that the smaller \$a{x, h 3 ) + 1|, the 
more accurate the approximation. 

These examples show that if one only needs a knowledge of the po- 
tential in the neighborhood of the oxide-semiconductor interface, the 
FSS provides a highly accurate approximation to the true solution 
much more cheaply than can be obtained by FD methods. In fact, to 
analyze three-phase devices, the cost of obtaining a FDS goes up 
sharply while the cost of a FSS remains nominal. For example, it took 
only 15 seconds and 34 K of memory to obtain the solutions presented 
graphically in Figs. 15 and 16. 

In reality, we are as much interested in the fields as we are in the 
potentials, and it is at this point that the difficulty with using the FD 
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method for solving these problems becomes most acute. In Figs. 6, 8, 
and 10, the dashed curves are plots of —(dtpm/dx) along y = hi, ob- 
tained from the FDS just discussed by differencing. In Fig. G for ex- 
ample, the two curves differ by nearly an order of magnitude at their 
peaks. If we take the FSS and difference it to estimate the first deriva- 
tive, we get a result which, in the neighborhood of the peaks, differs 
by at most 3 percent from the derivative obtained by differencing 
the FDS. We can conclude that the fields obtained from the FDS are 
badly in error. In order to calculate the fields from the FDS with any 
degree of accuracy, even for these simple examples, we would have to 
take a mesh so fine that the cost would become prohibitive. 
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